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We give a simple discussion of ghosts, unitarity violation, negative norm states and quantum vs 
classical behavior in the simplest model with four derivative action - the Pais-Uhlenbeck oscillator. 
We also point out that the normalizable "vacuum state" (in the sense defined below) of this model 
can be understood as spontaneous breaking of the emergent conformal symmetry. We provide an 
CO ■ example of an interacting system that couples the "particle" and "ghost" degrees of freedom and 

nevertheless remains unitary on both classical and quantum level. 

CN 

£ ' I. INTRODUCTION 

t CG ' 

The physics of systems with ghosts has recently attracted renewed attention The most interest in these systems 
is in connection with the theories of gravity. In particular the so called W 2 gravity, the theory with local conformal 
symmetry is known to have ghost modes |2| . This is usually considered to be a hindrance for a physical theory. Indeed 
an absence of a stable vacuum (lowest energy) state is disconcerting and is likely to lead to an instability, whereby 
the evolution extracts energy from the negative energy modes and pumps it into the positive energy modes producing 
a runaway instability. 

On the other hand, conformal gravity possesses much improved UV properties compared to the Einstein gravity, 
which render it renormalizable by power counting. The gravitational force in the theory of conformal gravity at large 
distances differs from the Newtonian gravitational force and this modification is capable of explaining data on galactic 
rotational curves with one fitting parameter, without introducing of the dark matter hypothesis (3j. It has also been 
suggested that conformal gravity may be able to solve the cosmological constant problem Q. 

These attractive features of the theory prompted attempts to solve the ghost problem. One approach attempts to 
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separate the ghost modes from the positive norm gravitons and ban their propagation "by hand" [1J . Another attempt 
is to quantize the theory using a nonstandard definition of a quantum mechanical normQ,0 following a more general 
program of quantizing PT invariant but non hermitian Hamiltoniansjy]. In a free limit this is essentially equivalent 
to treating the ghost modes as purely imaginary, which flips the sign of the ghost part of the Hamiltonian. It is not 
as yet clear whether any of these approaches can work in the full interacting theory. 

On the other hand, the instability in question may not be necessarily a fatal flaw. This is especially so in a theory 
of gravity, which governs the evolution of the universe and thus never actually relaxes to its ground state. Thus the 
nonexistence of a ground state in gravity may be just a way of life. In particular it has been suggested that a negative 
pressure due to ghosts may be a cause of the cosmological acceleration [8|. It has been argued that the time scale in 
which instability develops is way too short in theories which contain ghosts in the matter sector [jj. We are unaware 
however of a similar analysis of gravitational ghosts themselves, that is the ghost partners of the gravitons that arise 
in conformal gravity: the rapid decay of the vacuum discussed in|9( may be preventable if the ghost coupling to 
gravitons is nonlocal[l(|- It is not obvious therefore that the last word on viability of theories with ghosts has been 
uttered yet. 

The purpose of this note is rather modest and pedagogical. The potential interest notwithstanding, theories with 
ghosts are still considered somewhat esoteric and are not frequently discussed in particle physics literature. We aim to 
discuss pedagogically the simplest example of a theory with ghosts - the Pais-Uhlenbeck oscillator [ll|. Our goal is to 
explicitly demonstrate in this simple framework the meaning of some rather paradoxical notions that are sometimes 
used in the context of theories with ghosts, like negative norm states and violation of unitarity in theories with 
ostensibly perfectly hermitian Hamiltonian. We also demonstrate explicitly by solving the time dependent evolution 
in this theory how the soft UV behavior arises in dynamical context. 

We stress that the Pais-Uhlenbeck oscillator is in fact a unitary theory even though it possesses a ghost mode, and 
also give an example of a theory of interacting "particle" and "ghost" modes which is nevertheless unitary on the 
quantum level. All the above statements apply to quantum mechanical systems with the standard Dirac norm, as we 
do not recourse to a non standard quantization approach a la @. 

The Pais - Uhlenbeck oscillator was the subject of several papers in recent years, and its solution is well known^, 
12l [l3| . Nevertheless we feel that our simple and straightforward approach to the problem is illuminating and is worth 
recording. 
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II. THE PAIS-UHLENBECK OSCILLATOR 



The Pais-Uhlenbeck system is the theory of a single degree of freedom which satisfies a fourth order equation of 
motion. It is defined by the Lagrangian 

d 2 d 2 

For definiteness we assume u>\ > o>2- Our aim is to study the Hamiltonian dynamics with the view of quantum 
mechanical system, since a discussion of evolution of a wave function is most convenient in the Hamiltonian formalism. 
Although there exist a general formalism for calculating a Hamiltonian of four derivative systems, developed by 
Ostragradsky[14j, we find it more straightforward in the context of this particular model to introduce a pair of 
variables 

d 2 d 2 
X = —z + lo 2 iZ - Y = —z + lo 2 2 z (2.2) 

and consider them as independent coordinates. The rationale of this choice is, that the fourth order equation for the 
variable z 

(sM-° (2 - 3) 

can be written as a pair of second order equations for X and Y 

^X + ^X = 0; *Y + W ?Y = (2.4) 



To find the Hamiltonian we introduce the Lagrange multipliers for the constraints eq. (|2.2| 

d 2 d 2 
L = XY + a( — z + oj iZ -X)+ P(—z + u 2 z - Y) (2.5) 

Canonical momenta are calculated in the standard fashion pi = dL/dii. This definition leads to the following 
constraints 

P x = Py = 0; Va=V(i = -z; p z = -a - (2.6) 
The Hamiltonian, calculated in the standard way as the Legendre transform of the Lagrangian is 

H=~{p a +pp)p z -XY-a{ulz-X)-P{ulz-Y) (2.7) 

Commuting (calculating the Poisson bracket of) H with the primary constraints, eq. (|2.6|) we obtain secondary con- 
straints 

[H,P x ]=a-Y = 0; [H, Py] = /3 - X = 0; [H,p a - p p ] = (w? - u\)z - (X - Y) = (2.8) 
These can be used to express a, (3 and z in terms of X and Y, 

a = Y; = X\ z= X . - -= (2.9) 

The Dirac procedure for constraint systems requires that we use the Dirac brackets instead of the Poisson brackets to 
derive equations of motion. The net result of switching to the Dirac brackets is clear without a detailed calculation. 
The new "commutation relations" are such that the dynamical variables "commute" with all the constraints. Also, the 
modification is present only for those variables whose Poisson bracket with the original constraints does not vanish. 
Without any calculation the result in the present case is obvious 

Pa = K Y ; P/3=TTX; Pz = - ul){TTx - TTy) ( 2 -10) 

with the Dirac brackets 

[T H ,X j ] D = -6 ij (2.11) 
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The Hamiltonian then becomes 

H = -nA[n 2 Y - tt|] + -^Y 2 - -^X 2 (2.12) 
2 1 Y xi 20A 2flA v ' 

where we have defined Q = Ul +^ 2 ■ A = u)\ — W2- Finally rescaling the variables tt x = (fiA) 1 / 2 ^; x — (f2A) _1 / 2 X, 
and similarly for y we obtain 

= 2 71 * 2 lW ~ 2 7rx ~ 2 2 ' 
In terms of the new variables, the original coordinate z is expressed as 

1 



V4^A 



(x-y) (2.14) 



This is a very simple Hamiltonian. It is not bounded either from below nor from above, but nevertheless it generates 
a perfectly acceptable evolution. The two degrees of freedom x and y are decoupled, and classically each one simply 
satisfies a harmonic oscillator equation of motion. Classically there are no runaway solutions for these equations of 
motion for an arbitrary initial condition. Quantum mechanically the system also possesses finite positive norm states 
which evolve unitarily with time. 

Nevertheless the common jargon is that this quantum theory has negative norm states. In the next subsection we 
will clarify what this statement technically means, and stress that it is not a hindrance for a peaceful existence of a 
unitary evolution in this model. 



A. The Negative Face of a Divergent Integral 

The Hamiltonian of the Pais-Uhlenbeck system is not bounded from below. This is unusual and somewhat dis- 
turbing, since we normally expect that any open system will interact with some external degrees of freedom and 
generally relax to its ground state by loosing any excess energy to those degrees of freedom. However if the system is 
closed, no such loss of energy is possible and unboundedness of energy from below does not have to be a problem. In 
particular, in the present case the two harmonic oscillators do not interact with each other, no energy transfer from 
one to another occurs and the evolution is perfectly unitary, provided at the initial moment in time we start with a 
state which is localized at finite values of x and y. 

If one does insist, however to push the system to the lowest energy state, the evolution becomes non- unitary. This 
is simply due to the fact that this state is localized at infinite values of x, is non-normalizable and thus the probability 
"leaks" through the spatial boundary. 

To see this explicitly, let us define creation and annihilation operators in the standard way 

fu [Y 

a= \j2 X + l^ x (2 ' 15) 

The Fock vacuum of a is the normalized Gaussian state 

o|0) = 0; (x|0) = Ne-%* 2 (2.16) 

This is the state with highest energy in the x-sector. 

One can also formally define a state which corresponds to lowest energy eigenvalue, as the vacuum of a^[l2j 

a f |$) = 0; (x\$) = N_e% x2 (2.17) 

This state is non-normalizable and not physical, since a particle in this state is localized exclusively at infinity. The 
probability to find the particle at finite value of coordinate vanishes, since in the infinite volume limit the normalization 
constant N- vanishes faster than exponentially. 

Nevertheless in a certain formal way it corresponds to the lowest energy state. To see this, write the Hamiltonian 
for x mode in the standard form 

H x = -cjaat + E (2.18) 
Consider a tower of states above |$) generated by the action of operator a. 

H\$)=Eo\$); H\l) = Ha\<5>) = E a\<5>) - waa f a|$) = (E ... (2.19) 
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Thus applying operator a increases the energy of the state by ui, and the spectrum seems to be bounded from below. 
Another formal argument suggests that at least some of these states have negative norm. Let us calculate the norm 
of the "one particle state" 

(1|1) = ($|a f a|$) = ($|aa f - = -($|$) (2.20) 

Taken literally, this argument suggests that either the "one particle" state or the "vacuum" state has a negative norm. 
This is the origin of the usual statement that the theory has negative norm states. 

In fact, of course the norm of both of these states is positive once we regulate the system by putting it into a finite 
volume. The "vacuum" state is just a Gaussian which grows at large values of x. Its norm is positive in finite volume, 
and diverges (while remaining positive) as the infrared cutoff is removed. The one particle wave function can be found 
explicitly 

^ = {fr + {^i)^=^'' < 2 - 21 > 

The norm of this state obviously is also positive, and is even more divergent than that of the vacuum in large volume. 
None of the norms is negative. The flaw in the formal eq. (|2.21l) is of course precisely the fact that the states in 
question are not normalizable. To interpret the expectation value of a) a as the norm of a one particle state, one needs 
to act with a) on the bra, which amounts to integration by parts of the derivative in a) . The integration by parts 
however is not allowed, since the wave function grows at infinity. In particular 

<1|1)^H$)| 2 (2.22) 

as one can easily verify by an explicit calculation. In fact the difference between the two sides of the inequality is 
infinite. Thus "negative norm" is merely a jargon which refers to the fact that neither the norm nor matrix element of 
any reasonable operator like x n or p n is defined in the states of the form ea. (|2.19l) due to strong infrared divergence. 

Sometimes the procedure described above is referred to as a "quantization scheme", in the sense that the states 
of the tower eg. (12.191) do not belong to the Hilbert space of normalizable states. The unitarity in this quantization 
scheme is broken exactly for the reason explained above. All the wave functions with finite number of "excitations" 
above the "vacuum" |$) live on the edge of space. Once an infrared regulator (which makes the norm finite) is 
removed the wave functions vanish everywhere in the bulk. Such states run great risk of disappearing through the 
boundary under time evolution. 

On the other hand it is clear, that states which are created by the action of on |0) are normalizable and their 
evolution is perfectly unitary. One is normally interested in the situation when a particle can be detected in the bulk 
with finite probability. This physical condition makes the non-normalizable states physically irrelevant and devoid of 
interest. 



III. THE DEGENERATE CASE A = 

A special case of the Pais-Uhlenbeck system is when the two oscillators have the same frequency, A = 0. In terms 
of analogy with the W 2 gravity, this case is the most interesting. In this section we discuss some interesting features 
of the equal frequency limit. 

A. The Fate of the Normalized Wave Functions 

The limit A = of the previous expressions is a little tricky, since the transformation between the original variable z 
and x, y becomes singular. It is therefore not straightforward to take the limit directly on the level of the Hamiltonian. 
One cannot simply drop the terms in the Hamiltonian which naively vanish in the limit A — > 0, since the operators 
that multiply A may have divergent matrix elements. To illustrate this, let us first rewrite the Hamiltonian in terms 
of variables X and z, avoiding any singular redefinition of variables (here the variable x is defined as originally: 

X = £z + u, 2 z). 

H = —ttxtTz + X 2 - 2(fi + 7^A) 2 zX + 2Afl(n + ^A) 2 z 2 - ^Arhr| (3.1) 
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Suppose we naively drop the last two terms in ea. (l3.1|) . which formally vanish in the limit A — > 0. 

H = -^tt x tt z +X 2 -2(fl+^A) 2 zX (3.2) 

Let us now look for Gaussian eigenstates of the resulting Hamiltonian. Recall that at nonzero A we had four Gaussian 
eigenstates 



exp± l-^-I 2 ± -^-Y 2 \ = exp± I - X 2 ± 2cu 1 nAz 2 =F 2uj 1 Xz\ 

P \2flA 2flA J y \lo 1 ±lo 2 J 



(3.3) 



Three of these were non-normalizable and only one was the well behaved normalizable state peaked at x, z — 0: The 
normalizable state is 

* = exp-|^-X 2 + 2wif7Az 2 -2wiXz| (3.4) 

However if we seek all Gaussian eigenstates of the truncated Hamiltonian cq. (13.21) , we find only two states 

exp±|-^(X-2fi 2 z) 2 + 2f7 3 z 2 | (3.5) 

Evidently none of these two states is normalizable. These two Gaussian states are indeed obtained in the limit A — >• 
from two of the states eq. (|3.3p . Thus we seem to find no normalizable Gaussian eigenstates of a quadratic Hamiltonian 
cq. (|3.2j) . even though for any finite A a normalizable Gaussian eigenstate exists. This means that the Hamiltonian 
cq. p.2[) is not diagonalizable, which indeed can be formally proven [l^|.(7|. 

This conclusion is however a little hasty, as it is based on neglecting the last two terms in eq. (|3.1| . However, 
even though these terms are multiplied by A, in order to be able to neglect them, we need to be sure that they 
have vanishing matrix elements in the limit A — > 0. It is easy to see that this is not the case here. Indeed, in the 
normalizable state eq. (|3.4l) we have 



(z 2 )~(*x)~^ <■>■<»> 



so that in fact the last two terms in eq. Q3.ip are finite in the limit A — > and therefore cannot be simply discarded. 

The normalizable state eq. (|3.4[) does not disappear without a trace in the degenerate limit, but rather tends to a 
delta function of X 

^ 2 (X) -> S(X) (3.7) 

The action of the Hamiltonian eq. (|3.2[) on this state is ambiguous due to the first term in the Hamiltonian. One does 
obtain this state unambiguously, however as the equal frequency limit of eq. (I3.4[) 

Thus on the normalizable states the auxiliary variable X is frozen at zero, while the original variable z fluctuates 
freely with infinite amplitude. 

Interestingly, this suggests that in a sense the oscillator looses half of its degrees of freedom and also becomes 
"classical". Recall that the variable X is essentially the classical equation of motion for half the original modes of z, 
since X = ^z + uj 2 z. In the limit A — > 0, this quantity is fixed at zero without fluctuations. On the other hand the 
coordinate z itself fluctuates without restriction. Thus essentially the quantum system becomes a classical oscillator 
which can oscillate with arbitrary amplitude. 



B. Dynamical conformal symmetry 

As an interesting aside, we note that at A = the theory dynamically develops a conformal symmetry, which 
is spontaneously broken by normalizable states. For the purpose of this discussion, it is convenient to revert to 
normalization in which the Hamiltonian is simplest in the limit A — > 0, eq. (|2.13|) . In the equal frequency limit the 
Hamiltonian cq. (|2.13|) is invariant under the following transformation 

x — > a; cosh t + t/sinhi; y — > ycosht + xsinhi; z — > e~*z (3-8) 

It is natural to refer to this symmetry as conformal. This symmetry is not obviously present in the Lagrangian 
cq. ([2.1[) . In fact the Lagrangian is multiplied by a constant under the transformation eq. (|3.8p . However, as we have 



6 



seen in the previous subsection, in the equal frequency limit the dynamics of z is such that on normalizable states it is 
pinned to satisfy X = -^z + uj 2 z = 0. As a result the Lagrangian vanishes for all physically interesting configurations. 
Scaling of the Lagrangian by a finite factor therefore is indeed a "dynamical" symmetry in this limit. 

Interestingly this symmetry is spontaneously broken, in the sense that the normalizable "vacuum" , or in fact 
any of the normalizable physical states, is not invariant under it. The wave function of the "lowest energy" , the 
non-normalizable eigenstate of the operator is indeed invariant under the conformal transformation: 




However for the normalizable Gaussian 




(3.10) 



It is clear that any state whose wave function is localized at finite values of x and y is necessarily not invariant 
under the transformation eq. (|3.8[) . Thus the conformal symmetry is "spontaneously broken" on normalizable states. 
Since the representations of conformal group eq. (|3.8p are infinitely dimensional, the finite energy spectrum is infinitely 
degenerate. This is of course well known and obvious since adding any number of excitations of the x oscillator and 
the same number of excitations of the y oscillator docs not change the energy in the degenerate limit (l^. It is 
nevertheless amusing, that this degeneracy can be understood as a spontaneous breaking of conformal symmetry. 



IV. DYNAMICS: CLASSICAL VS QUANTUM 



The dynamics of the classical Pais-Uhlenbeck oscillator is identical to that of two decoupled harmonic oscillators. 
The variables x and y satisfy the harmonic oscillator equations of motion, and the fact that the energy of the x- 
oscillator is negative is irrelevant, since the energies of each oscillator are separately conserved. 

Quantum mechanically, however the situation is very different. Here the overall sign of energy is reflected in the sign 
of the phase of the wave function. For the evolution of states which are initially product wave functions 1 ^i(x)^>2(y) 
this is again unimportant, however it affects strongly the time evolution of "entangled" states. The simplest calculation 
where the quantum mechanical importance of the sign flip for the x-oscillator manifests itself, is the propagator of the 
z. It is of course well known, that the UV behavior of the propagator in four derivative theories is much softer than 
in theories with ordinary kinetic term. The Pais-Uhlenbeck oscillator is the simplest example of this kind. Although 
this is a trivial calculation, we present it here for completeness. 

A. The propagator 

To calculate the propagator of z we need to calculate the propagator of x and y separately. For y this is the usual 
harmonic oscillator calculation. 

1. The y propagator 

The Hamiltonian for the y mode is 

H=±p 2 + ~u iy 2 (4.1) 

The annihilation operator a 

a = a I 

2 va uii 
evolves in time according to 

a(t) = a(0)e-' lUlt (4.3) 



a=M/£(y+^) (4-2) 
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For the Feynman propagator: 



G y (t) =< T{y(t)y(0)} > = — < Q(t)[a(t)a* (0) + a\t)a(0)} + 9(-i)[a t (0)o(t) + 0(0)0* (t)] > (4.4) 



wc have 



G (t) = [G^e^ 1 * + e(-t)e ZWlt ] (4.5) 

To perform the Fourier transform, as usual we introduce the regulator which makes the integral converent for large 
times 



G y (p) = — / rfie lpt [9(i)e-^ 1 *e- £ * + 9(-t)e laJlt e £t ] = — TT — (4.6) 



1 

This is the standard result, which upon integration over the frequency p gives the equal time expectation value in the 
vacuum 

<' 2 > = /f G » = i (47) 

2. The x propagator 

The propagator of x is equally easy to calculate in the physically relevant "vacuum"- the highest energy state. The 
Hamiltonian now is 

H=-\p 2 -\u 2 x 2 (4.8) 



and 



a=J — (£+—); oft) = a(0)e lW2t (4.9) 

2 w 2 



The same calculation as before now gives 



G x (t) = (0|T{a:(t)a;(0)}|0) = -L[e(t)e iW2 * + 6(-t)e-^*] (4.10) 

2uj 2 



and 



G x {p) = ~ 2 — (4.11) 

p z — U>2 — It 

This differs from cq. (l4.6[) by the overall sign and also by the sign of the regulator e. As is easily seen, these two sign 
changes cancel each other in the calculation of equal time quantities. For example 

(0|* 2 |0) = / Z Gx{p) = ^ ( 4 - 12 ) 

which is the correct result for the normalizable Gaussian eigenstate of the x oscillator. 

3. The z propagator 

Finally combining the results for x and y, and noting that due to the symmetries of the system the mixed propagator 
vanishes (x(t)y(0)) — 0, we obtain 

Again, this is the standard result, showing a softened UV behavior, since the propagator of z vanishes much faster 
for high frequencies than that of a harmonic oscillator. This indicates of course, that the time evolution of z is very 
smooth and has a very small high frequency component. 
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4- The "propagator" in the unbounded state 

What happens if we try to calculate the propagator of the x oscillator in the unbounded Gaussian state? Of course, 
as explained above this calculation is purely formal, since the integrals over this wave function are divergent. Still, 
formally proceeding as before we can define 

G-{t) = {*\T{x(t)x(Q)}\$) (4.14) 
We still use eq. (|4.9[) . but this time it is cv that annihilates the state <£■. We then formally obtain: 

G~(t) = -7^-[Q(t)e- lul2t + e(-t)e ZUJ2f } (4.15) 

and 

G-(jp) = 2 ~* . (4.16) 

The sign of the regulator e is now the same as for the positive energy harmonic oscillator, which is simply the reflection 
of the fact that the state |$) is formally the lowest energy state of the system. However this propagator leads to the 
same paradox of negative norm states as discussed in the previous section. Calculating the equal time expectation 
value, which should be by definition positive, we find 

(^ 2 |<I>H-^ (4-17) 

This again underscores the point, that non-normalizable states, if manipulated formally, can be mistaken to have 
negative norm. 

B. Time evolution: the wave function 

It is instructive to see explicitly how the wave function of the system evolves in time. In particular we would like 
to see the origin of the smooth UV behavior of the Pais-Uhlenbeck system in terms of the time evolution of wave 
functions. 

We are mostly interested in the degenerate case A — > 0, and will therefore study time evolution generated by the 
Hamiltonian 

1 d 2 1 d 2 1 9 , 1 , , , 

2 dy 2 2 dx 2 2 y 2 K ' 

We want to follow the time dependence of simple quantum averages, like (z 2 (t)) and ((X(t) + Y(t)) 2 ). The first 
observable is the obvious choice, since it is the fluctuation of the coordinate of the original oscillator, while the second 
one is the fluctuation of the second order equation of motion. We will choose an initial state such that both these 
operators have sensible (finite) averages. 

We are not interested in states which are simple product states of the form ipi(x)tp2(y)- As far as the expectation 
values of all Hermitian operators go, the evolution of such a product state is identical to that of a state tpi(x)ip2(y) 
evolved with the positive energy harmonic oscillator. We will thus be interested in states which are not trivial product 
states in the variables x and y. A simple initial wave function that satisfies these requirements is 



r(0) = \>x,»^-i 



Afl , N2 1 t N2 

.An i . 2 .Af2 i , 2 .An i . 
{ -e + 4n^A )x +( ^ + 4n^A )y +2( ^~4nVA )xy 



(4.19) 



Note that we have scaled out the dependence on the frequency difference A explicitly. Strictly speaking for non- 
vanishing A we also have to keep the frequencies of the two oscillators in the Hamiltonian different. However the 
Hamiltonian itself is smooth in the degenerate limit, and it is only the relation between x, y and z that involves 
divergent coefficients. Thus with the appropriate choice of the wave function we can make z finite also at A — > 0. 
Specifically, for the state ea. (|4.20p we have 



(z 2 )=r 2 : 



({X + Y) 2 )=£ 2 (4.20) 



9 



Since the evolution is free, a Gaussian wave function preserves its Gaussian shape at any later time. Thus at any time 
t we have 



-~A(t)x 2 - ~B(t)y 2 - C(t)xy 



ip(t) = N(t) exp 

Acting on this wave function with the Hamiltonian we obtain the evolution of the coefficients 



(4.21) 



A = i\A 2 - C 2 - n 2 ]; B = i[C 2 — B 2 — fi 2 ]; C = iC[A-B). 
After some algebra this leads to 



(4.22) 



C = C 



A + B 
A + B 



which is solved by 



with 



C{t) = a[A{t) + B(t)} 

C(0) 
A(0)+B(0) 

Using this result for C(t) in eci. (|422|) . and defining + B(t) = u{t), A(t) - B{t) = v(t). we have: 

u = iuv 



i>^[(i-2a> 2 + i« 2 -2ft 2 ] 



with the initial conditions: 



t(0) = 2 



= /AO 



1 



V £ 2 40r 2 A 



ft is easy to see that the solution has the form 



u(t) = 



1 



/+ + /_ cos 2flt 



, «(*) = 



w(0) = 



-£2/_nsm2fii 
/+ + /_ cos 20i 



(4.23) 



(4.24) 



(4.25) 



(4.26) 



(4.27) 



(4.28) 



(4.29) 



where f± are constants determined by the equations of motion and the initial conditions. After some algebra, for the 
initial conditions eq. (|4.28p we obtain 



and 



A(t) = 
B(t) = 
C{t) = 



/ ± = §( ± i + ^V) g2+4 ^ A2r2 



SI (e 2 + 4Q 2 A 2 r 2 ) - z2A(l - Q 2 £ 2 t 2 ) sm2SU 
2A (1 - cos 2flt) + n 2 £ 2 T 2 (l + cos 2Qt) 

n (£ 2 + 40 2 A 2 t 2 ) + »2A(1 - n 2 £ 2 T 2 ) sin 2Qt 
2A (1 - cos 2flt) + n 2 £ 2 T 2 (l + cos 2Qt) 
n e~ 40 2 A 2 r 2 



(4.30) 



(4.31) 



2A (1 - cos 2Qt) + 2 £ 2 r 2 (1 + cos 2fit) 
The time dependent probability density can be written as: 



i/j'ijj = N 2 exp 



Q 2 40 2 A 2 r 2 n 

'2A (X + V ' (1 - cos 2VLt) + fl 2 £ 2 T 2 (l + cos 2Qt) ~ 2A^ X ~ V ' 



e 



(1 - cos20i) + 2 £ 2 r 2 (l + cos 20t) 

(4.32) 
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Thus we find 



< z 2 (t) > = - 
((X(t)+Y(t)) 2 ) = 1 



1 



n 2 e 
1 

ft 2 r 2 



(1 - cos2Qi) + t 2 (1 + cos2fii) 
(1 - cos 2ttt) + £ 2 (1 + cos 2f2t) 



(4.33) 



These expressions are notable for their absence of features. Normally one expects that if the initial state is very far 
from the vacuum, the evolution should delocalize it in a short time, so that the amplitude of the fluctuation of the 
coordinates should become very large. This is exactly what happens in the standard positive Hamiltonian harmonic 
oscillator, as we will demonstrate in the next subsection. However eg. (14.331) shows that in the Pais-Uhlenbeck system 
both interesting averages evolve smoothly in time on the scale determined by the initial state averages. Clearly, if 
both t and £ are finite, the averages stay finite throughout the evolution. This is despite the fact, that the "vacuum" 
of the system is such that r 2 oc 1/A— s-oo, £ 2 cx A — >• 0, as discussed in the previous section. If we start the system 
"close" to its vacuum state, that is with £ 2 cx 1/t 2 cx A, it is still true that at all times parametrically the averages are 
the same, fluctuation with the amplitude proportional to the initial average. Thus it does not matter, if the system 
starts off far from the vacuum, or close to it, the evolution is smooth and the averages at all times are proportional 
to those in the initial state. 

To underscore that this is very different from the standard harmonic oscillator, we perform the same exercise as 
above for the two decoupled oscillator systems. 



C. The baseline: oscillators with positive energy 

We now consider time evolution generated by 

1 f> 2 1 Ft 2 1 1 

H = --^~-^ + -n 2 i/ + -n 2 x 2 (4.34) 

2 dy 2 2 dx 2 2 J 2 V ' 

For a Gaussian wave function eq. (|4.21jl the evolution of the parameters A(t), B(t) and C(t) is given by: 

A = i[-C 2 - A 2 + n 2 ]; B = i[-C 2 - B 2 + fi 2 ]; C = -iC(A + B) (4.35) 
This is simplified for our initial state where A(t) and B{t) stay equal for all times: 

A = i[-C 2 - A 2 +n 2 }; C=-2iAC, (4.36) 



with initial conditions given as in I4.20I This is solved by 



c= 1 A= igtosm&m + t) 

f + g cos (2Clt + 4>Y f + g cos (2Clt + <j>) ' ' 



provided 

1 

Imposing the initial conditions, 

AQ 1 igflsind) „. . AQ I I 

A(O) ^ + 4^A^ 2(/ + g cos0) ; C(0) = ^-4^A^7T^- (439) 

we find 



f - g 2 = 7^2 ( 4 - 38 ) 



£ 2A n 2 T 2 e~i . , ^ i4r» 2 A 2 r 2 +e 2 , ^ nA n 2 T 2 e + i , AAn , 

where $ = i<f>. Finally, the solution for our initial conditions is 

n (40 2 A 2 t 2 + £ 2 ) cos2M + i2A(n 2 T 2 £ 2 + l)sin2ftt 



A(t) = Bit) 



2A Q2 T 2 £ 2 (1 + cos 2fit) - (1 - cos 2Qt) + i(2n 2 Ar 2 + J^) sin2Qi 



o 4Q 2 a 2 t 2 - f 2 

C(t) = 77^ ^-=-^ a (4.41) 

2A n2 T 2^2(! + cos 2 Qt) - (1 - cos 2fit) + i(2ft 2 Ar 2 + ^) sin 20* 
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For small A we expand these expressions to second nontrivial order 



A(t) = B(t) 



C(t) 



iQtan2Q,t + 



sin 2Qt 



2QA (n 2 T 2 £ 2 - 1) cos 2VLt + (n 2 r 2 £ 2 + 1) cos Ant 

1? sin 2 2ttt 

2f2A ft 2 r 2 £ 2 (l + cos 2Q£) - (1 - cos 2Qt) 

~fF sin 2 2ttt 



(4.42) 



Generically at arbitrary time we have 



Re[A + C] oc Re[A-C] cx 



A 



(4.43) 



sin 2 2flt 



and thus 



((x-y) 2 ) cx ((x + y) 2 ) cx 



sin 2 2flt 



(4.44) 



A 



This is precisely what one normally expects. Our initial state is very far away from the ground state. It was chosen in 
such a way that the center of mass coordinate x + y had large fluctuations, 0(1/ A), whereas the relative coordinate 
x — y had small fluctuations 0(A). One expects a state like this to expand very quickly and become delocalized in 
all coordinates. Indeed eg. (14.441) displays precisely this feature: the relative coordinate fluctuates with amplitude of 
order 1/ A almost all the time, except for a very short time interval St cx A within every period of evolution. 

Thus indeed, we see that the time evolution of the Pais - Uhlenbeck oscillator is smoother than that of a system of 
decoupled harmonic oscillators, in the sense that the averages in the Pais- Uhlenbeck case fluctuate on the scale given 
by the initial state and do not develop additional large variations throughout the evolution. 



We have seen that the quantum evolution of the Pais-Uhlenbeck oscillator is unitary. This is not very surprising, 
nor very exciting since the two second order degrees of freedom in this case are decoupled, and each one follows a 
Harmonic oscillator evolution. In fact the system has two conserved quantum numbers - not just the total energy, 
but also the energy of each individual oscillator is conserved. For this reason the classical motion in the X, Y plane 
is bounded and the quantum evolution is unitary. 

A more interesting and general question is whether interacting systems with ghosts can be unitary. The worry is 
clear. We have a Hamiltonian which is unbounded neither from above nor from below, and once the two modes x and 
y are allowed to interact, there is a real and present danger that the system can develop an instability, where both x 
and y run away to infinity even though the total energy stays conserved. 

In the quantum mechanical context one can pose the following question: does a system of coupled "particle" and 
"ghost" degrees of freedom possess normalizable eigenstates. If the answer is affirmative, such system enjoys unitary 
quantum evolution, since the probability to find the system in finite volume does not decrease with time [l 71 ] . If this 
is not the case, such systems would not allow for unitary quantum mechanical evolution and probability would leak 
out completely through the boundaries in a finite amount of time. 

The aim of this section is to present a simple example of a model, which remains unitary even though it contains 



For definiteness we choose /i > 0. At /i = the theory is clearly unitary, as the particle and ghost degrees of freedom 
are decoupled, and evolution of each one separately is unitary in exactly the same sense as for the Pais - Uhlenbeck 
oscillator. 

The question about stability can be asked already on the classical level. It was noted in [16| and also [13J, that 
some systems of this kind allow for classically stable solutions, namely oscillatory solutions for which the amplitude 
does not grow without bound as a function of time. Specifically ref. [la ] studied numerically the evolution of eq. (|5.ip 
for X12 = and found that the classical behavior of the system is stable as long as the initial energy stored in the 
oscillators is not too large. Denoting the initial displacement of the oscillators from the equilibrium by M, ref.flU] 
found that for M 2 < M 2 — j^ui 2 the behavior is oscillatory, while for M 2 > M 2 the amplitude of oscillations grows 
without bound. The addition of the quartic self interaction Ai j2 further stabilizes the system. We have repeated the 
numerical exercise of [l6j for the system eq. (|5.1|) . and have found a similar behavior in a wider range of parameters. 



V. 



A SIMPLE UNITARY INTERACTION 
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In fact as long as the coupling fi remains small fj, -C Ai,2 we did not see classical instability for any initial conditions 
that we have tried. Examples of evolution for several initial conditions are given in Fig.l. This suggests that when 
the interaction is weak enough, the classical system is absolutely stable, although it is not possible to prove such a 
statement by numerical methods. 




(c)cco = 1, yo = 1, xo = 0, yo = {d)xo = 1, yo = 1, xo = 0, yo = 



FIG. 1: Typical time evolution of x (red) and y (blue) for different initial conditions. The parameters are chosen as cji = 
3, Ai = 10, u)2 — 5, A2 = 7, fi = 3. The evolution is plotted over two time intervals to show the detailed structure of time 
dependence and to demonstrate the absense of instability over very long times. 

Note, that in order for the quantum system to be unitary, its classical counterpart has to have stable evolution for 
arbitrary initial conditions. Otherwise quantum tunneling will connect stable and unstable regions of the phase space 
and will inevitably lead to violation of unitarity. This is the situation, for example in the upside down Mexican hat 
potential U(x) = —X(x 2 — x^) 2 . Classical solutions with total energy — Xxq < E < and initial displacement \x\ < xo 
are regular. However quantum mechanically the system is non-unitary due to finite probability of tunneling into the 
unbounded region \x\ > xq. 

In the present case, one can give an argument that the theory remains stable, at least in a limited range of param- 
eters. Let us consider the limit uj\ < W2. In this case one can use the classical Born-Oppenheimer approximation. 
Since y oscillates much faster than x, one can consider the motion of y in the background of fixed x. Thus for given 
x the dynamics of y is given simply by an anharmonic oscillator with the frequency, which for large x behaves as 
u! 2 = \ix 2 . This is clearly a well defined bounded motion. The dynamics of x is affected by the average value of y 2 
for a given trajectory. Given the initial energy E stored in the mode y, we have (for large x, which is the interesting 
and potentially dangerous regime) y 2 oc Ej \ix 2 . The dynamics of x then is governed by the effective potential 

\<jj\x 2 + X 2 x 4 - [iy 2 x 2 = ^uj 2 x 2 + A 2 x 4 - E (5.2) 

Thus the dynamics of x in this approximation is unaffected by y and is bounded and stable. A similar argument can 
be given for the opposite case uji 3> w 2 ■ Thus at least when the two frequencies are very different there is no instability 
for arbitrary initial conditions. In this case one expects that the quantum theory is well defined and unitary in the 
sense explained above. 

In the next subsection we present another line of reasoning supporting the same conclusion for small /i. 
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A. Asymptotics of Eigenfunctions for Small fi 



One way to establish that a quantum theory has normalizable eigenstate is to find asymptotics of eigenfunctions 
for large values of coordinates x and y. 
As usual, we introduce an eikonal S via 



- M P - s ( x <y) 



* = Ne 



(5.3) 



If the eikonal is positive and divergent for large values of the coordinates, the wave function is normalizable. For large 
values of S, \x\ and \y\ it satisfies the following "semiclassical" equation: 



2 [dy~ 



1 fdS_ 

2 { dx 



) 2 + A l2/ 4 



2„,2 



X 2 x + ijlx y 



(5.4) 



We will not attempt to solve this equation in full generality, but rather explore the behavior of S for small values of 
ji. For jjL = the solution is simply a sum of the solutions for two decoupled degrees of freedom: 



So(x,y) 



/2X 



2 i i3 

x\ d 



(5.5) 



The crucial point is that the structure of the potential is such that for fi <C A,;, the perturbation is smaller than the 
leading order potential for generic large values of x and y. This is of course very different from the standard pertur- 
bation theory around a harmonic oscillator potential, where a perturbation is usually bigger than the unperturbed 
potential for large values of the coordinate. Thus although the standard perturbation theory around a Harmonic 
potential is asymptotic, we expect the perturbation theory in /i to have a finite radius of convergence. 

Let us therefore solve eq. (|5.3p perturbatively. Let S = Sq + Si, where Si oc fi. We first solve the equation for 
x, y > 0. To first order in fi we have: 



Changing variables x 



l 



/2\ 2 x 



, y 



V2A7y 2 ^ + V^x 2 §± + mxV = 
ay ox 

and defining x — x ± y the equations becomes simple 



dSi 4/i 



dx~ Xi\ 2 (x+ 2 - x- 2 ) 2 



A well behaved solution to this equation is: 



(5, 



(5.7) 



Si(x + , x 



4/, 



1 



AiA 2 2x+ 3 



-2 _ rf+2 



arctanh 



(5.8) 



In terms of the original variables, the solution can be written as: 



S 1 (x>0,y>0) = V2^ 2 - 2 ^y-^ 2X 



ifix y 



(VATy + v%x) 2 
Extending the solution to other regions of the plane we find 



x 3 y 3 



— logU; ) 



Si(x,y) = V2^x 2 y 2 



\%\y\ - v%\x\ 
(VX~i\y\ + Vx-2\x\) 2 



kl 3 bi 3 locr ATM, 

(VAT|y| + V^\x\) 3 &K \ \2\x\> 



(5.9) 



(5.10) 



As expected, the correction Si is smaller than So at large values of the arguments, and thus the asymptotics of the 
wave function is determined by So- Thus we find that for small /i our model quantum mechanically has normalizable 
eigenstates, and therefore unitary evolution. 
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